Helical scattering and valleytronics in bilayer graphene 
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We describe an angularly asymmetric interface-scattering mechanism which allows to spatially 
separate the electrons in the two low-energy valleys of bilayer graphene. The effect occurs at 
electrostatically defined interfaces separating regions of different pseudospin polarization, and is 
associated with the helical winding of the pseudospin vector across the interface, which breaks the 
reflection symmetry in each valley. Electrons are transmitted with a preferred direction of up to 60° 
over a large energetic range in one of the valleys, and down to —60° in the other. In a Y-junction 
geometry, this can be used to create and detect valley polarization. 

PACS numbers: 72.80.Vp, 73.22.Pr, 85.75.-d 



I. INTRODUCTION 

Research on low-dimensional systems has received a 
significant stimulus by the recent isolation of atomically 
thin sheets of graphite, known as graphene. 1 Graphene's 
attractive electronic properties can be traced back to 
a unique low-energy band structure (corresponding to 
massless Dirac particles in pristine monolayers, which be- 
come massive in bilayers), whose key features are dictated 
by the existence of an intrinsic orbital degree of freedom, 
known as pseudospin — a choice of two n orbitals, pre- 
sented by the two chemically identical carbon atoms in 
the unit cell of the graphene honeycomb latticed The 
celebrated phenomenon of chirality ties the pseudospin 
to the momentum, which results in a Berry phase of it 
in monolayers (2ir in bilayers) whose consequences are 
observed in the quantum Hall effector— and in the chiral 
Klein tunneling across np junctions.^' 7 

A descended, remarkably robust electronic feature — 
shared by graphene monolayers, bilayers, and multilay- 
ers of various stacking order — is the splitting of the low- 
energy band structure into two energetically degenerate 
but orbitally inequivalent valleys (of different isospin), 
situated around the corners of the Brillouin zone (the 
K and K' points). In this work we demonstrate that 
the electrons belonging to the valleys can be separated 
and detected by a variant of the Klein tunneling set-up, 
which utilizes a certain type of electrostatically defined 
interfaces in bilayer graphene. The interfaces in question 
separate regions of opposite pseudospin polarization, and 
are obtained by combining pairs of oppositely charged 
top- and back gates as illustrated in Fig. Ufa). In the 
past, such interfaces have attracted attention because 
they guide valley-dependent states^ that propagate along 
the interface; for transport in the perpendicular direc- 
tion, they realize a pseudospin variant of a spin valve^i 
which allows electronic confinement Ai By resolving the 
transport across the interface angularly in each valley, we 
find that the direction of perfect transmission is skewed 
away from normal incidence, see Fig. HJc)-(h). The bro- 
ken reflection symmetry is intimately related to a hidden 
helical structure — the winding of the pseudospin vector 
as one crosses the interface, which is opposite in the two 



valleys as illustrated in Fig. |TJb). The ensuing helical 
scattering favors transmission at opposite angles of inci- 
dence, and provides a means to separate electrons in the 
valleys by directing them to different leads. 

From a practical perspective, efforts to create an imbal- 
ance in the population of the valleys are motivated by the 
prospects of electronic applications that do not directly 
rely on the charge of the carriers. 12 Dubbed valleytronics, 
this scheme aims to implement analogies to the rapidly 
maturing field of spintronics,— with the two valleys sub- 
stituting the two states of physical electronic spin. The 
original proposal to create and detect valley polarization 
relies on edge effects j^ 2 - alternative proposals employ bro- 
ken sublattice symmetry combined with in-plane electric 
fields^ or drastic doping^ Besides the practical obsta- 
cles to implement these routes, these proposals face the 
severe problem of strong intervalley scattering occuring 
at pn junctions in monolayer graphene. 1 ^ In this context, 
the valley-dependent helical scattering offers two key ad- 
vantages: The required interfaces can be implemented 
and controlled by ordinary gates, and intervalley scat- 
tering is negligible. Towards the end of this paper, we 
evaluate the feasibility of this route to valleytronics by 
numerical results for wide (but finite-sized) Y-junction 
devices. 

The paper is organized as follows: Section [TTI presents 
analytical results for scattering at pseudospin interfaces 
in bilayer graphene. Section Mil presents numerical re- 
sults for valleytronics applications. Section IIVI contains 
conclusions. Appendices provide background for the an- 
alytical and numerical methods. 



II. HELICAL SCATTERING 

We start our considerations with the scattering at 
pseudospin interfaces in bulk graphene, which allows a 
fully analytical treatment that reveals the connection of 
valley-dependent scattering and helical winding. 

In bilayer graphene, pseudospin refers to the probabil- 
ity amplitude to find an electron on the top layer (spin 
up, amplitude V , t( r )) or on the bottom layer (spin down, 
amplitude ipi(r)). Pseudospin polarization can be ere- 
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FIG. 1: (Color online). Helical scattering of electrons in bi- 
layer graphene at an interface between regions of different 
pseudospin polarization. Panel (a) sketches in a side-view 
how such an interface is obtained by sandwiching a bilayer 
graphene flake between positively and negatively charged 
gates (orange and blue, respectively; not drawn to scale). 
Panel (b) shows the rotation of the pseudospin vector in the 
valleys around the K and K' points. In the polar plots of pan- 
els (c)-(h), solid (blue) lines show the angle-resolved transmis- 
sion probability T(tp) in the K valley, while dashed (red) lines 
show this probability in the K' valley. Each panel pertains to a 
different Fermi energy Ef, ranging from values just above the 
gap where the transmission in each valley is angularly asym- 
metric but small [panels (c) and (d)], over a regime where 
angular asymmetric scattering occurs around a transparent 
direction [panels (e)-(g)] to the regime far above the gap, 
where one approaches valley-independent, angularly symmet- 
ric scattering [panel (h)]. Panels (i) and (j) contrast these 
finding to the case of Klein tunneling at an np junction, 6 ' 7 
which does not distinguish the valleys. 



ated by oppositely charged top- and back gates, which 
induce an electrostatic potential [7-j- = — t/j, = U(x) in 
the two layers. Patterning of the gates allows to spa- 
tially vary the potential. The relevant low-energy elec- 
tronic transport in this landscape can be described by a 
two-component Hamiltonian 17 



H 



U(x) 
-(fix +i£p v ) 2 /1rn 
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(1) 



which acts on the spinor wave function ip(r) — 
(-0t( r ) I '04.( r )) T - Here in is the effective mass, and fi x , 
fi v are the two components of the momentum operator in 
the graphene plane. The valley index £ takes the values 
£ = 1 in the K valley, and £ = — 1 in the K' valley. 

A sharp pseudospin interface along the line x = is de- 
scribed by a step-like profile of the electrostatic potential, 
with U(x) = —Uo for x < and U(x) = Uq for x > 0. 
In the regions of constant potential, the electronic bi- 
layer spectrum is gapped, so that electrons can only move 
freely for Fermi energies Ep fulfilling \Ep\ > A/2 = \Uq\. 
(At energies within the gap, there is a pair of states which 
propagate along the interface^ these states do not con- 
tribute to the transport across the interface. In the fol- 
lowing, we will assume Ep > Uq > 0.) In each valley, the 
angle-resolved probability T(ip) = \t(ip)\ 2 of electrons ap- 
proaching the interface from the left to be transmitted to 
the right can be calculated in a standard wave-matching 
procedure (see Appendix [TffJ. The transmission ampli- 
tude is 
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" — 2m(E F -Uo)' H ~ 2m(E F -U 

Uq ) 1 / 4 . The resulting angle-resolved probability only de- 
pends on Ep/Uo, and is plotted in Figure [TJc)-(h) as a 
sequence of polar plots. 

Just above the gap (Ep — 1.02 Uo), the transmis- 
sion probability in both valleys is small for any angle 
of incidence (this is the essence of the pseudospin valve 
effect^). As the Fermi energy increases, the transmission 
probability increases, but in an angularly asymmetrical 
fashion: over a large range of energies, electrons in the 
K valley have large transmission probability for incident 
angles in the range of 30° to 60°, where the interface 
is almost transparent. The transmission characteristics 
in the K' valley are similarly skewed, but into the op- 
posite direction. An approximately symmetric, valley- 
independent transmission probability only emerges for 
energies of the order of Ep — 2C/ and beyond. These 
observations can be quantified analytically by consider- 
ing the direction of unit transmission probability, which 
exists for values Ep > ^/9/8f/o ~ 1.06 Uq, and then is 




given by the momentum vector 



Po 



( 3 ) 
The component parallel to the interface is finite, and has 

an opposite sign in the two valleys; it only becomes neg- 
ligible when Ep ^> [To- 
rt is natural to inquire whether there is a fundamen- 
tal mechanism behind this valley-dependent transmission 
pattern. Here we offer two observations — firstly, the he- 
lical winding of the pseudospin vector for propagating 
modes as one crosses the interface, which manifests the 
broken reflection symmetry in each valley; secondly, a 
condition for the pseudospin vector of evanescent modes, 
which identifies the transparent direction p . 

The pseudospin vector P = (tr) is a unit vector com- 
posed of the expectation values of the Pauli spin matri- 
ces. For a constant potential U, a propagating wave with 
momentum p carries pseudospin 
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The out-of-plane component is a measure of the pseu- 
dospin polarization, while the projection onto the 
graphene plane is linked to the propagation direction 
(this is the celebrated chirality of charge carriers in 
graphene 2 ); notably, the sign of the y component de- 
pends on the valley index £. Far to the left and to the 
right of the interface, the pseudospin vector has a finite 
but opposite out-of-plane component. As U changes from 
— Uq to Uq the pseudospin vector describes an arc on the 
unit sphere, which has opposite orientation in the two 
valleys (see Fig. 0Jb)), and thereby explicitly breaks the 
reflection symmetry in a fixed valley. The asymmetric 
transmission pattern is therefore connected to a hidden 
helical structure in each of the two valleys. 

Our second observation involves the evanescent modes, 
which carry pseudospin close to the interface, and whose 
amplitudes are related to the propagating modes by the 
wave matching conditions. At the transparency condi- 
tion ([3]), the evanescent modes on both sides of the inter- 
face have the same amplitude and phase on both layers, 
■)/>cv oc (1, 1) T . Their pseudospin vector P cv = i therefore 
lies in the graphene plane, and points into the direction 
perpendicular to the interface. While this condition is 
the same in both valleys, the properties of the evanes- 
cent modes differ; as a result, they meet the condition 
for the finite, valley-dependent value of p y = po y speci- 
fied in Eq. ([3]) , which again explicitly breaks the reflection 
symmetry in each valley. 



III. VALLEYTRONICS APPLICATIONS 

We now turn to the prospects to employ helical scat- 
tering for valleytronics applications. In view of previous 
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FIG. 2: (Color online). Absence of intervalley scattering. 
Panel (a): Top- view sketch of a (zigzag) graphene nanorib- 
bon geometry with a pseudospin interface, and definitions 
of reflection and transmission coefficients (summed over all 
channels) when a K-polarized current is injected into lead 
f. Panels (b) and (c): Numerical results for the valley- 
resolved reflection and transmission coefficients (nanoribbon 
width W ~ 150h/y/2m\U \). 



proposals in monolayer graphene ) 12 ' 14 i 15 a key question 
is whether pseudospin interfaces induce intervalley scat- 
tering in zigzag nanoribbons, as this is the orientation in 
which the valleys are conserved away from the interface. 
This question goes beyond the continuum description, 16 
and is numerically addressed in Figure [2] The top panel 
shows the nanoribbon geometry. To quantify the inter- 
valley scattering, consider the current l[^ in lead I = 1,2 
and valley X=K, K' in response to a bias voltage V 
applied to lead I, where the electrons are polarized in 
valley K (below we describe how such valley polariza- 
tion can be created; lead 2 is grounded). For small 
bias we can define valley- and lead-resolved conductances 
G 1 ^ = I l {x/V, which can be evaluated using standard 
quantum transport algorithms (see Appendix [Bjl . The 
bottom panels show valley-resolved transmission coeffi- 
cients Tx = G\x/Gq (where Go is the conductance quan- 
tum) , as well as reflection coefficients Rx , defined analo- 
gously based on the reflected electron stream, for a wide 
nanoribbon of width W ~ lhQ%/ yj2m\Uo\ (for these wide 
ribbons we do not find any difference between the zigzag 
and antizigzag edge configurations^). The results show 
that intervalley scattering is negligible both in reflection 
as well as in transmission. Encouragingly, the transmis- 
sion coefficient exceeds the reflection coefficient once the 
regime of helical scattering is entered. The results also 
reveal that for the assumed sharp interface the reflec- 
tion coefficient is almost energy-independent (softening 
the interface further reduces the reflection and interval- 
ley scattering). 

In a suitable geometry, the helical scattering should 
therefore enable to spatially separate the electrons in the 
two valleys. To evaluate the feasibility of such a bulk val- 
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FIG. 3: (Color online). Valleytronics applications. Panel (a): 
Top-view sketch of a valley splitter, consisting of a Y-junction 
where electrons are injected through lead 1, and collected in 
leads 2 and 3. In this set-up, the contact region has the 
same gate polarity as lead 1, and the incoming current in 
this lead is valley- unpolarized. Panels (b) and (c): Valley- 
resolved conductance of electrons in leads 2 and 3. Panel (d): 
Ratio of the larger to the smaller conductance in these two 
leads. Analogously, panels (e)-(h) demonstrate that valley 
polarization can be created in lead 1 by injecting unpolarized 
currents into lead 2 or 3. Panels (i)-(l) show how the total 
current in leads 2 and 3 can be used the detect the valley 
polarization of an incoming current in lead 1. 



ley filter we consider a Y-junction geometry with three 
leads, shown in Figure[3fa), where a pseudospin interface 
is placed at the entrance of leads 2 and 3. This is done 
with the intent that electrons injected into lead 1 are 
preferably transmitted into lead 2 when they are in the 
K valley, and into lead 3 when they are in the K' valley. 
To quantify the anticipated valley polarization, we now 
consider the current l\ x in lead I = 2, 3 and valley X=K, 
K' in response to a bias voltage V applied to lead 1 (the 
injected current from this lead is now valley-unpolarized; 
leads 2 and 3 are both grounded), and the associated 
conductances G lx = l[ x /V. Figures EJb)-(d) show the 
energy dependence of these conductances, along with the 
ratio of the larger over the smaller conductance in each 
lead. The results demonstrate that a finite valley polar- 
ization can indeed be obtained for conditions where the 
conductance itself is also large (Ep ~ I.IUq, where the 
preferred transmission directions are « ±60°). By inject- 
ing unpolarized currents into leads 2 or 3, the Y-junction 
geometry can also be used to create valley polarization in 
lead 1; this is shown in panels (e)-(h). Panels (i)-(l) show 
how valley polarization in lead 1 is detected by the total 
current in leads 2 and 3, which is quantified by conduc- 
tances G\ x (now X stands for the valley polarization of 
the incoming current). A K valley-polarized current in- 
jected into lead 1 will preferably divert to lead 2 (panel j), 
while a K' valley-polarized current will preferably divert 



FIG. 4: Alternative design of the valley splitters, polarizers 
and detectors of Fig. El where the contact region has the same 
gate polarity as leads 2 and 3. Qualitatively, the results for 
valleys K and K' are now interchanged; however, the influence 
of edge states between leads 2 and 3 result in a notable lead- 
asymmetry of the conductance ratios. 



to lead 3 (panel k). 

Figure |4] shows an alternative design of the Y-junction, 
where now the contact region has the same gate polarity 
as leads 2 and 3. This produces larger conductance ratios 
for valley polarization and detection, but suffers from 
larger asymmetries which originate from edge states that 
connect leads 2 and 3. Symmetry upon interchanging the 
lead index is restored when one simultaneously inverts 
the sign of U and E F . 



IV. CONCLUSIONS 

In summary, we demonstrated that electrostatically 
defined pseudospin interfaces in bilayer graphene result 
in angularly asymmetric, valley-dependent scattering. 
Above a certain energy threshold, the transmission be- 
comes perfect into a particular direction, which is oppo- 
site in the two valleys. This direction only turns towards 
normal to the interface when energies are very large. The 
broken reflection symmetry in each valley can be made 
manifest by considering the helical winding of the pseu- 
dospin vector across the interface. 

From a conceptual point of view, our findings enforce 
the role of interfaces to unravel fundamental properties of 
charge carriers in graphene, thereby complementing the 
earlier instances of np junctions (which allow to realize 
Klein tunneling^— and Veselago lensingi£), flake bound- 
aries (which can support edge states^), and regions of 
large charge carrier density (which induce surface states 
that support tunneling over large distances^). From a 
more practical perspective, the helical scattering opens 
an avenue to manipulate the valley degree of freedom in 
bulk graphene. 



The author gratefully acknowledges discussions with 
Elsa Prada, Pablo San- Jose, and Edward McCann. 



Appendix A: Wave matching 

The transport problem across the pseudospin interface 
can be solved by matching the propagating and evanes- 
cent modes to the left and right of the interface, which 
are obtained from the Schrodinger equation Hip = Eip 
with Hamiltonian H given in Eq. ((1|). 

The scattering at the interface conserves energy E and 
the component p y of the momentum parallel to the inter- 
face. In the region x < 0, where U — —Uq, an incident 
mode with fixed p y and E has longitudinal momentum 



p x = J2my/E 2 — U^ — p 2 = p, and is described by a 
wave function 



ip in = cxp(ipx + ip y y) 
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FIG. 5: (a) Tight-binding model of a flake of bilayer graphene, 
consisting of a top and a bottom layer where each layer ac- 
commodates a honeycomb lattice of carbon atoms. Black dots 
denote dimer sites, where carbon atoms of both layers lie on 
top of each other and are coupled with inter-layer coupling 
constant 71. Light (red and blue) dots denote sites of iso- 
lated carbon atoms (either on the top or the bottom layer). 
Lines denote intra-layer nearest-neighbor hopping with cou- 
pling constant 70. (b) Effective model for the two low-energy 
bands, obtained after decimation (algebraic elimination) of 
the dimer sites. Thick lines denote coupling with strength 
27', while thin lines denote coupling with strength 7', where 

i = 70/71 ■ 



where a 



(P+iJPv) 
2m(E-Uo) 



The reflected wave 



V> rc f = cxp(-ipx + ip y y) i a „ j (A2) 

has longitudinal momentum p x — —p. There is also an 
evanescent mode 

V>icft = exp(qx + ip y y) J, (A3) 

which decays (to the left, away from the inter- 
face) with imaginary momentum p x = —iq, q = 

yj2my/W=U* + $, and /3 = ^^Ty 

In the region x > 0, where U = Uq, the transmitted 
wave 



^trans = exp(ipx + ip y y) ( 1 



a" 



(A4) 



has longitudinal momentum p x = p, and the decaying 
evanescent mode is given by p x = iq, 



bright = exp(-qx + ip y y) ( ° 



(A5) 



The wave matching conditions at the interface x = 
between the regions are continuity, 



iWiMij-'ffMfj.w 



and continuity of the first derivative, 



'"I 1 a)- pr (a*)- iqc ( I) =pt {l )+^ d (i 



Together, they form an inhomogeneous system of four 
linear equations which uniquely determines the four co- 
efficients r, t, c, and d, where in particular 



t = 



(A7) 



2ip 2 y {l - /? 2 )Ima + 2ipq[Rea(l + (3 2 ) - /3(1 + |a| 2 )] 
p2(l - P 2 ){1 - a* 2 ) + ipq[(a* - f3) 2 + (1 - a*p) 2 ] ' 

The transmission probability follows from T = \t\ 2 (the 
corresponding reflection probability is-R=|r| 2 = I — 
T). (In the main text, we renamed E — > Ep, and p — > 
p x .) Condition ([3]), which singles out the propagation 
direction po with unit transmission T = 1 , coincides with 
the condition /3 = 1. The transmission amplitude then 
takes the value t = (a — !)/(! — a*). 



Appendix B: Numerical Method 

We base our numerical calculations on a tight-binding 
model which describes the two low-energy bands of bi- 
layer graphene. This model is obtained from the physical 
bilayer graphene lattice of carbon atoms, shown in Fig. 
[S{a), where the 7r orbitals of neighboring atoms in the 
same layer are coupled by an intra-layer hopping energy 
70 , while orbitals of atoms that lie on top of each other 
(dimer sites) are coupled by the inter- layer hopping en- 
ergy 71 . At low energies, the latter atoms are only visited 
in virtual transitions, and can be eliminated using deci- 
mation. This amounts to solving the Schrodinger equa- 
tion for the amplitudes on the dimer sites, and inserting 
the result back into the remaining equations of the non- 
dimer sites. The (generally energy-dependent) self en- 
ergies and mediated hoppings can be set to their values 
at the Dirac point E = 0. The resulting two-component 
tight-binding model is shown in Fig.[SJb). Nearest neigh- 
bors are now coupled with strength 27', where 7' = 70/71 



is an effective hopping energy. Coupling to next-nearest 
neighbors is still constrained to atoms of different sub- 
lattices, and carries strength 7'. This model allows to 
efficiently enter the regime of many propagating modes 
in the leads without interfering with the split bands. 

The valley-resolved conductances are determined from 
the relation 



1 — <J"0 



E 

elcad /, valley X 
melcad 1 



the factor 2 arises from the degeneracy of the physical 
spin, and S nm are components of the scattering matrix, 
which we calculate by a standard numerical decimation 
method^ The conductances Gj x and G L 1X are defined 
and obtained analogously. The valleys are resolved by us- 
ing the Bloch wave propagation factor A = cxp(iv3afc x ) 
of modes in leads aligned along the valley-preserving 
zigzag orientation (here a is the intra-layer distance of 
neighboring carbon atoms). In the K valley, ImA < 0, 
while in the K' valley ImA > 0. 



Here Go = 2e 2 /h is the conductance quantum, where 
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